Critical Behavior of the Two-Dimensional Randomly Driven Lattice Gas 
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We investigate the critical behavior of the two-dimensional randomly driven lattice gas, in which 
particles are driven along one of the lattice axes by an infinite external field with randomly changing 
sign. A finite-size scaling (FSS) analysis provides novel evidences that this model is not in the same 
universality class as the driven lattice gas with a constant drive (DLG), contrarily to what has been 
recently reported in the literature. Indeed, the FSS functions of transverse observables (i.e., related 
to order-parameter fluctuations with wave vector perpendicular to the direction of the field) differ 
from the mean-field behavior predicted and observed within the DLG universality class. At variance 
with the DLG case, FSS is attained on lattices with fixed aspect ratio and anisotropy exponent equal 
to f and the transverse Binder cumulant does not vanish at the critical point. 
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Phase transitions are characterized by a drastic change 
in the macroscopic behavior of many-body interacting 
systems when control parameters are varied. In the case 
of critical phenomena the onset of the ordered phase is 
accompanied by fluctuations on all length scales. In spite 
of the difficulties in accounting efficiently for such cou- 
pled fluctuations, we have currently a deep and detailed 
understanding of the collective behavior of a large class of 
systems in thermal equilibrium thanks to a variety of re- 
sults and methods, both analytical and numerical. Crit- 
ical collective behaviors, on the other hand, are also ob- 
served in the steady states of systems far from thermal 
equilibrium ]J. In contrast to equilibrium cases, the sta- 
tionary probability distribution is not known a priori and 
the possible occurrence and nature of a phase transition 
can no longer be determined by usual entropy-energy ar- 
guments. The absence of detailed balance, the generic 
algebraic decay of space-dependent correlations as func- 
tions of the distance, their strongly anisotropic scaling 
properties, and the flux of particles or energy through the 
system are among the general features which make these 
problems particularly difficult and rich in phenomenol- 
ogy. Because of the lack of a general framework, it is still 
worth focusing on specific toy models in order to gain in- 
sight which might possibly lead to a more comprehensive 
theory. Among the models characterized by nonequilib- 
rium steady states the simplest is the uniformly driven 
lattice gas 2] (DLG), a generalization of the Ising model 
to nonequilibrium conditions due to the action of an ex- 
ternal nonconservative force, inducing a particle current 
through the system. Although the DLG was introduced 
more than 20 years ago, there is still room for debate on 
the nature and properties (in particular the universality 



class) of the nonequilibrium critical behavior observed 
upon changing the temperature. At first, the relevant 
feature of the model was assumed to be the presence of 
a particle current However, more recently, this point 
has been criticized by arguing that the strong anisotropy 
is, instead, its qualifying, character 0,0- In addition to 
the theoretical debate |6J, seemingly contradictory evi- 
dences are also coming from Monte Carlo (MC) simu- 
lations 0, S H, 13 • Tne I sm S model can be driven to 
nonequilibrium conditions also by the coup ling to two 
thermal baths at different temperatures [lfj, controlling 
the hopping rates of the particles in different lattice di- 
rections. In a simpler version, the temperature affect- 
ing the jumps in one direction is taken to be infinite. 
Hereafter we will refer to this direction as the longitudi- 
nal one (||) whereas to the remaining as transverse ones 
(_L). This model is equivalent to a DLG in which the 
microscopic external driving force is along the longitu- 
dinal direction, with infinite modulus and a sign that 
is randomly chosen for each lattice site every time step 
(annealed randomness). The resulting model is called 
the randomly driven lattice gas (RDLG) [nj. At vari- 
ance with the DLG, no net particle current is flowing 
through the system. MC simulations |l2| indicate that 
in the RDLG transverse fluctuations of the order pa- 
rameter (i.e., the fluctuations with wave vector in the 
transverse direction) are not effectively described by a 
Gaussian model. Indeed, in two dimensions, the case 
we shall consider from now on, MC simulations give 
u± = 0.62(3) and n = 0.13(4) 12]. These estimates 
rely on a field-theoretical result for the anisotropy expo- 
nent A = (z'u/V-l) - 1, i-e., A = 1 - r]/2 wl Q which 
enters the finite-size scaling (FSS) Ansatze used to ex- 
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tract critical exponents. In Ref. |f| the numerical data 
for the RDLG and the DLG on the same finite lattices 
have been compared. According to these data the two 
models have the same finite-size scaling (FSS) behavior. 
If true, this implies that they belong to the same univer- 
sality class and thus the strong anisotropy and not the 
particle current is the relevant feature in determining the 
leading critical behavior of driven diffusive systems. The 
same conclusions has been drawn in Ref. Q by studying 
the short-time dynamics, although some points of the 
analysis therein remain unclear |13( . Here we reconsider 
the problem presenting the results of a new series of MC 
simulations of the RDLG. The critical behavior of the 
system (proper of the thermodynamic limit) is extracted 
from data on finite lattices by means of a FSS analysis 
that does not require free parameters ^J, in contrast 
with that previously employed . We briefly recall the 
definition of the RDLG. Consider a rectangular lattice 
and for each site x introduce an occupation variable n Xl 
which can be either zero (empty site) or one (occupied 
site). The external field E is acting along the longitu- 
dinal direction but with random sign. The dynamics of 
the model is of Kawasaki type: A lattice link (xy) is 
randomly chosen, and, if n x ^ n y , a particle jump is 
proposed and then accepted with Metropolis probability 
w(PAH + (3Ei), where £ = (1, 0, —1) for jumps (along, 
transverse, opposite) to E, w(x) = min(l,e _;r ), and 
AH is the variation of the standard lattice-gas nearest- 
neighbor attractive interaction H = —4 J2{ xy ) n x n y due 
to the proposed jump. The parameter (3 plays the role 
of an inverse temperature. In the DLG, E is constant 
and time-independent. Periodic boundary conditions in 
the direction of E make it a nonconservative field and 
drive the system into a nonequilibrium stationary state. 
Although the boundary conditions are not so relevant in 
the RDLG, we will assume them periodic in all directions. 

At half filling, the RDLG undergoes a second-order 
phase transition. Indeed, at high temperatures the 
steady state is disordered whereas at low temperatures 
the system orders: The particles condense forming a strip 
with interfaces parallel to E. These two phases are sep- 
arated by a phase transition occurring at the critical 
value (3 C (E) depending on the field E. Here we will 
concentrate on the particular case in which E is infi- 
nite. For a strongly anisotropic system in d dimensions, 
with finite size Lii x L^ 1 , the FSS limit corresponds to 
t = 1 — (3/f3 c — > (where (3 C is the bulk critical tem- 
perature), Li\,L± — > co, keeping fixed both combina- 
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and therefore also the so-called 



aspect ratio 5 a = Ll/^ 1+ ^ /L± [J^j. Accordingly, the 
FSS analysis of numerical data generally requires an a 
priori knowledge of the exponent A [17| . It would be 
a real step towards a better understanding of nonequi- 
librium phase transitions to have at disposal FSS in a 
form suitable for these systems, reliable and powerful 



enough to disentangle those key features which might be 
buried in tiny differences when the volume of the sam- 
ples are increased. In this direction we have already 
performed a detailed study of the FSS of the DLG M 
by using the general strategy introduced in Ref. |14| . 
confirming the mean-field behavior of transverse fluctu- 
ations, with A = 2, in agreement with the predictions 
of Ref. 0- It is therefore a crucial test to examine by 
the same method also the RDLG. For previous studies 
of FSS in strongly anisotropic systems see Refs. I3E3- 
We will show that the method is sensible enough to high- 
light the differences in the critical behavior of the DLG 
and RDLG (contrarily to the claims in Ref. to an 

extent that goes beyond the numerical differences in the 
critical exponents and probes the spatial structure of cor- 
relations. In spite of the generic power-law decay of the 
two-point correlation function (n x no) for large x [12j, it 
is possible to define a finite-volume correlation length 
0. Given the Fourier transform G{q) of (n x no), one 
considers G±(q) = G({<7|| = 0, q± = q}) and defines the 
correlation length 
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where q n = 2sin<7„/2 is the lattice momentum and 
q n = 2nn/L±. Hereafter we will denote L± simply by 
L. In Fig.^we report the FSS plot of the ratio £,2l/£,l, 
where &l and £l are computed at the same temperature 
but lattice sizes 2L and L, respectively, keeping constant 
the aspect ratio with A = 2, i.e., 5 2 ^ 0.200 0. For 
comparison we report as a solid line the mean-field pre- 
diction which is approached by the DLG data on the 
larger lattices 0- In the present case, deviations from 
the mean-field behavior increase with increasing lattice 
sizes. Note that, if 5i is the correct aspect ratio for the 
model, then one observes the crossover towards the FSS 
of the model in the strip geometry L± = oo, when keep- 
ing S2 constant and fixed [ljj. Accordingly, the 
points corresponding to larger lattices in Fig. ^ eventu- 
ally accumulate on some limiting curve as L increases, in 
agreement with the predictions based on the field- 
theoretical model of Ref. 11]. Fig. □ refers to geometries 
with A = 1, i.e., Si « 0.223 (upper set of points) and 
Si ~ 0.326 (lower set). Note that we used A = 1, al- 
though a correction of the order of rj/2 to this value is ex- 
pected. In these cases we have been able to reach L = 88. 
In contrast with the mean-field behavior, which does not 
depend on the specific value of 5a , now we do expect a 
dependence of the FSS curves on the actual value of Si. 
The critical properties can be extracted from the previ- 
ous plot as follows. For a given 5i, consider the scaling 
function £,2l/£,l = F(z) as a functions of z = £l/L and 
expand it around z* , which is defined as the point such 
that F(z*) — 2, i.e., as the value of £l/L at the critical 
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FIG. 1: FSS plot of the transverse correlation length £l in 
the two-dimensional RDLG with fixed S2 ~ 0.200. Crosses, 
circles, squares, triangles, diamonds, and stars correspond to 
lattices of increasing size L = 14, 16, 18, 20, 22, 24. The Gaus- 
sian behavior approached by the DLG is given by the solid 
line. 
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FIG. 2: FSS plot of the transverse correlation length £l in the 
two-dimensional RDLG with fixed Si = 0.223 (upper set of 
points) and Si = 0.326 (lower set). Empty squares, empty tri- 
angles, diamonds, crosses, full circles, full squares, full upright 
triangles, and full downright triangles correspond to lattices 
of increasing size L = 20, 22, 24, 28, 32, 36, 40, 44. 



temperature. Denoting Sz — z — z*, one finds 9] 
F{z) = F(z*) + F'(z*)Sz + 0[{5zf] 

= 2 + ^(2 1 /^ -l^)Sz + 0[(5z) 2 } . 

A linear fit of our data gives z* — 0.1337(3) for S\ — 
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FIG. 3: FSS plot of the ratio £,l/xl in the two-dimensional 
RDLG with fixed Si = 0.223 (upper set of points) and 0.326 
(lower set). Symbols are as in Fig. 01 



0.223 and z* = 0.1594(1) for Si = 0.326, with the same 
u±_ = 0.61(3). The corresponding critical temperature is 
the same in the two cases. Note that z* for St_= 0.326 is 
almost equal to the mean-field value l/(2n) |f|. Indeed, 
we had chosen this value of Si in order to be very close 
to the mean-field predictions and test whether the FSS 
method employed is able to detect the differences. These 
results suggest that at variance with the DLG, where 
mean-field scaling at fixed S2 is observed, in the present 
case scaling is attained only at fixed Si and is not com- 
patible with mean-field behavior. Indeed, not only does 
z* depend on the geometry, but also the critical exponent 
v±_ differs from the Gaussian value 1/2. The qualitative 
dependence of z* on Si is accounted for [2(J by the field- 
theoretical model of Ref. ^l|. A similar analysis, with 
similar results, has been performed for the susceptibility 
Xl = G±(qi) [23|. We find more instructive to present 
data for the ratio Al = Ci/Xii which is independent 
of £l/L (f° r L large enough) whenever the critical ex- 
ponent 77 vanishes. In Fig. [3] we present the FSS data 
of this observable, for the two values of Si previously 
considered. At variance with the DLG, where we got 
A2l/Al ~ 1, here we see a pronounced and systematic 
dependence on £l/L and on the actual value of Si- For 
the two different values of z*(Si) we do find the same 
value for AilIAl ~ 1.05, which is equal to 2 V [£| and 
leads to the estimate 77 = 0.07(1). Further evidence of 
differences in the critical behavior of the DLG and of the 
RDLG is provided by the transverse Binder cumulant 
g L = -G^\qi,qi,-qi,-qi)/[N L Gj_{qi)], where is 
the Fourier transform of the four-point connected correla- 
tion function at a given (3 and lattice geometry, computed 
at the first allowed transverse momentum q±, and Nl is 
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FIG. 4: FSS plot of the transverse Binder cumulant gh in the 
two-dimensional RDLG with fixed Si = 0.223 (upper set of 
points) and Si = 0.326 (lower set). Symbols are as in Fig. 



the total number of spins in the lattice. The FSS plot in 
Pig.|3|shows clearly that for both the values of Si consid- 
ered, g2h/9L — 1 at the corresponding z*{S\). Therefore, 
gL at the critical point does not vanish in the thermody- 
namic limit, at variance with what has been found for the 
DLG |9(. The estimated values of the critical exponents 
are in good agreement with previous numerical findings 
(althou gh our result for rj is smaller than that reported 
in Ref. [12() and theoretical estimates based on the field- 
theoretical model of Ref. Moreover, the universal 
FSS functions for the correlation length agree with those 
computed in field theory at first order in an e expansion 
around the upper critical dimension d = 3 |2(ij |. In con- 
clusion, we have shown that the FSS approach as devised 
in Ref. is sensible enough to distinguish clearly be- 
tween the critical behavior of the DLG and of the RDLG, 
two systems which in particular geometries may exhibit 
quite similar behavior for relatively small volumes and 
not too close to the critical temperature. Therefore, on 
one side we have a sound numerical method to exam- 
ine also nonequilibrium critical phenomena, on the other 
we have eventually established that the key features of 
these two models are different, as they do not belong 
to the same universality class. Th e ag reement with the 
field-theoretical analysis of Refs. |3J, [llj, [l2| suggests that 
indeed the leading critical behavior of the DLG is gov- 
erned by the presence of a particle current whereas that 
one of the RDLG is dominated by strong anisotropy. 
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